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ABSTRACT 

IC 10 X-l is a bright (L x =10 38 ergs -1 ) variable X-ray source in the local group 
starburst galaxy IC 10. The most plausible optical counterpart is a luminous 
Wolf-Rayet star, making IC 10 X-l a rare example of a Wolf-Rayet X-ray binary. 
In this paper, we report on the detection of an X-ray orbital period for IC 10 X-l 
of 34.4 hours. This result, combined with a re-examination of optical spectra, 
allow us to determine a mass function for the system f(M) = 7.8M and a 
probable mass for the compact object of 24 — 33M Q . If this analysis is correct, 
the compact object is the most massive known stellar black black hole. We 
further show that the observed period is inconsistent with Roche lobe overflow, 
suggesting that the binary is detached and the black hole is accreting the wind of 
the Wolf-Rayet star. The observed mass loss rate of [MAC92] 17- A is sufficient 
to power the X-ray luminosity of IC 10 X-l. 



Subject headings: stars: Wolf-Rayet — galaxies: starburst — X-rays: binaries - 
X-rays: galaxies 
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Introduction 



Models of the evolution of High Mass X-ray Binaries (HMXB) pred ict the existence 
of helium star+compact object pairs (Ivan den Heuvel fc de Loord Il973l ). Such systems 
should form at the very end of the X-ray binary evolution, after the secondary (donor) 
star has been stripped of it's hydrogen either through Roche lobe overflow or mass loss via 
a strong wind. In either ca se, a luminous helium star with a compact companion is formed 
( Ergma fc Yungelsonl I1998T ). These systems are expected to be rare. lErgma fc Yungelson 
( 119981 ) predict ~100 helium star+black hole pairs in the Galaxy. However, only a small 
number of these are expected to form accretion disks and hence be visible as X-ray sources. 
Identification of such rare systems is important, because they have the potential to put 
strong constraints on the evolution o f massive binary pairs . There are currently only three 
candidates: Cyg X-3, NGC 300 X-l JCarpano et alil2007ah and IC 10 X-l. 



IC 10 X-l is a bright ( L.,.=10 38 ergs l ) variable X-ray source in the local group metal- 



poor starburst galaxy IC 10 (IBrandt et al. 



1997 



Bauer &: Brandtll20041) . It is surr ounded by a 



shell of non thermal radio emission jYang fc Skillmanlll993h and X-ray emission (jWang et al.l 



20051 ; IBrandt et al.lll997l ) which may be associated with the supernova which produced the 
compact object in IC 10 X-l. There are 4 possible optical counterparts to t he X-ray source, 



with the most plausible being a bright Wolf Rayet star [MAC92] 17- A (jCrowther et al. 



20031 ). Spectroscopic observations of [MA C92] 17- A reveal promin ent He II line emission, 
suggesting an identification as a WNE star (IClark fc Crowtherll2004l ). In this paper we report 
the discovery of an X- ray orbital period in IC 10 X-l using data from SWIFT and Chandra 
( IPrestwich et al.ll2006l ). In Section [2] we describe the Chandra and S WIFT observations, and 
in Section [3] we discuss the constraints the period puts on the accretion mechanisms and 
mass of the compact companion. 



2. Obervations and Data Analysis 

IC 10 X-l was observed with Chandra on 2006 2 November and 2006 4 November 
for approximately 45 ks per observation (ObsIDs 7802 and 8458). Data were processed 
using standard CI AO analysis software, version 3.4. After extraction of a lightcurve, we 
observed a large (factor of 7) flux increase in the first Chandra observation followed by a 
similar flux decrease in the second observation (see Figure [T|) . The large flux change and 
sharp profile of the flux modulations suggested we were observing eclipses in IC 10 X-l. 
We immediately applied for, and were generously awarded, a SWIFT target-of-opportunity 
(TOO) observation of 100 ks. IC 10 X-l was observed with SWIFT beginning 2006 21 
November for approximately 700 seconds per 90 minute SWIFT orbit for a total of 97 
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ks spanning 246 hours. Data were processed by the SWIFT pipeline and screened using 
standard procedures outlined in the XRT Data Reduction Guide. 



2.1. Timing Analysis 



In order to search for periodic signals in the SWIFT observations of IC 10 X-l, we 
merge the SWIFT event files using the CIAO tool dmmerge. We extract the events in a 60" 
radius around the Chandra source position; this is reasonable due to both the extent of the 
S WIFT point-spread function and the sparsity of nearby sources as observed in the Chandra 
observations. The A Lomb-Scargle (LS) periodogram of IC 10 X-l is shown in Figure [2J The 
periodogram shows no evidence for increased noise at low frequencies ("red noise"), indicating 
that the LS method can be used to search for peaks in the power spectrum. The full data 
range was searched from the Nyquist frequency (twice the 2.5073 s XRT readout rate) to 
half the full observation duratio n of 123 hours. The nor malization used in the LS method 
is the total variance of the data ( Home fc Baliunasl 19861). Th e LS p eriodogram determines 
a period of 34.40. Following the method of lliorne fc Baliunas ( 1986 ) we determine a period 
uncertainty of ± 0.83 hours. 

Two sets of simulations to were carried out to ascertain the false alarm probability 
(FAP) of the signal. The first was a Monte Carlo simulation, in which fake spectra were 
generated with Gaussian noise set to the standard deviation of the actual data. The second 
was a bootstrap method, which creates a fake periodogram by randomly rearranging real 
data values. Observation times were taken from the real data, and 1 x 10 5 fake spectra 
were generated for each method. The number of simulated spectra with peak amplitudes 
greater than the the observed peak amplitude (i.e. the FAP) was in both Monte Carlo and 
Bootstrap simulations. Thus the FAP < 1 x 10~ 5 , and the significance > 4.5a 

IC 10 X-l was observerd in 2003 by both Chandra and XMM-Newton. The 30 ks Chandra 
observatio n of IC 10 X-l showed t he source to be variable, but did not show any evidence for 
an eclipse (IBauer fc Brandtll2004l ). A sharp increase in the X-ray flux o f IC 10 X-l durin g 
the 45 ks XMM-Newton observation is probably due to an eclipse egress (jWang et al.ll2005l ). 
The data in these observations are consistent with the period derived in this paper. 



3. Discussion 



The most plausible explanation for the regular flux modulations is an orbital eclipse of 
the X-ray emitting object by the donor star. It is unlikely the observed period is super-orbital 
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e.g. due to a precessing warped accr etion disk) because it i s shorter than most super-orbital 



periods (most are 10-100 days, e.g. Idarkson et al.l (120041 ) ). While we cannot rule out an 



alternative explanation without a fully sampled optical radial velocity curve, we assume the 
modulation is orbital for the rest of the paper. 



3.1. Constraints on the Accretion Mechanism 



The period can be used to constrain the accretion mechanism, in particular to determine 
whether the donor fills its Roche lobe, or whether accretion occurs via a wind. Here we 
test the hypothesis that the accretion is via Roche lobe overflow by determing whether the 
observed period is consistent with theoretical values of the period for compact object masses 
in the range 1-100M Q . 

The orbital period, P, is related to the mass of the primary (Mi), the mass of the donor 
(M 2 ) and the separation between them (a) via Kepler's Third Law: 



P 



2na 3 / 2 
[G(M 1 + M 2 )]V2 



Here we assume that the inclination of the system is close to 90 deg, as implied by the 
existence of an eclipse. If we assume that the donor fills its Roche lobe, we can make the 
simplifying assumption that the Roche lo be radius (Rjx) is equal to the donor radius {Rq). 
The separation can then be expressed as (lEggletonlll983l ): 



D 0.6g 2 / 3 + ln(l + g 1 / 3 ) , oN 
a = R2 0A9^ (2) 

where q = MijM\. If the mass and radius of the donor star are known, equations CD and 
|2]can be used to determine values of the period for plausible masses of the compact object. 



The mass of [MAC92] 17- A derived from spectroscopic data is ~ 35 M (IClark fc Crowther 



20041 ) . The uncertainty on this value is large. It is possible (but unlikely) that [MAC92] 17- A 
might have a mass as low as 17 M Q . The rad ius of [MAC92 ] 17- A is derived using standard 
Mass- Radius values for WR stars taken from Langer ( 19891 ). Assuming that [MAC92] 17- A 



is indeed the donor star, we find that for M 1 =1-100M Q the period derived using equations 
[Hand [2] is 2-3.5 hours. These conclusions still apply if the mass of [MAC92] 17- A is closer 
to 17 M Q than 35 M Q . Clearly, the observed period is inconsistent with Roche lobe overflow 
and implies that accretion is via a wind as seen in many other high mass X-ray binaries. 
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3.2. The Mass of the Compact Object 

The standard met hod to measure the mass o f star s in a binary system is to determine 



the mass function (e.g. iMcClintock fc Remillardl (120061 )): 



fiM ) = ™l- (3) 
f[M) - 2nG ~ (1 + g) 2 [6) 

where K 2 is the half amplitude of the velocity curve of the secondary. The mass function 
is the minimum mass of the compact primary. In the case of the IC 10 X-l system, we have 
reasonable estimates for the the orbital period and the inclination but require a value of 
K 2 to secure the mass function. The centroid of the Hell 4686A line is commonly used to 
determine the radial velocity of WR+OB binaries because it is formed in the inner wind, 
close to the star. In the light of the disco very of the orbital perio d in IC 10 X-l, we re- 



examined the optical spectra presented by IClark fc Crowtherl (120041 ) to search for shifts in 



the centroid of the Hell 4686A line which might be due to the orbital motion. We note that 
the spectrum of [MAC92] 17- A is characteristic of a WN star, and does not show evidence for 
features seen in accretion driven outflows (e.g. HI, Hel, Fell). We are therefore reasonably 
confident that the Hell 4686A line in [MAC92] 17- A is associated with the donor star and 
has minimal contribution from any accretion disk. 

Individual Gemini GMOS spectra of [MAC92] 17- A were re-extracted relative to the 
IC 10 nebular [OIII] 5006.8 A line. Figure [3] shows the individual spectra together with 
Gaussian fits to the Hell 4686A line. Two one-hour exposures were taken sequentially on 
the night of December 22 2001, and two further spectra obtained on the nights of January 
13 2002 and January 15 2002. As might expected for a binary system with a period of 34.4 
hours, there is no evidence for a shift in the velocity centroid of the stellar Hell 4686A line 
within IC 10 for the two exposures taken on December 22 2001. However, the line centroid 
on 15 Jan 2002 is definitely shifted substantially relative to the other datasets. The FWHM 
of the fitted gaussian for the 15 Jan 2002 dataset is narrower than that determined for the 
other epochs (10A vs 15-17A). This may be due to a CCD defect which contaminates the 
red component of the 15 Jan 2002 line. Our best estimate for the centroid of the Hell 4686A 
line on 15 Jan 2002 is ~ —300 km s -1 , compared to ~ +450 km s" 1 for the other epochs. 

We make the assumption that the shift in the Hell 4686A line centroid of approximately 
750 km s _1 is due to the radial velocity of [MAC92] 17- A. We further assume that 750 km 
s _1 corresponds to the maximum possible velocity displacement. We find that K 2 = 375 km 
s _1 and the mass function f(M) = 7.8 M . This value of the mass function implies that the 
compact object is a black hole. 



We can further constrain the mass of the compact object by using equation [3] and 
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assuming a range of reasonable values for the mass of the donor s tar and inclination of the 
system, as shown in Table HJ The mass of [MAC92] 17- A derived by lClark fe Crowtherl (120041 ) 
from spectroscopic data is ~ 35 M . There is considerable uncertainty in this estimate, and 
so we use a conservative lower limit of 17 M . The existance of deep eclipses suggests that 
the inclination of the system is close to 90 deg, however we cannot rule out the possibility 
that the inclination is considerably less and use 45 deg as a lower limit. Tabled] demonstrates 
that the minimum mass of the black hole to satisfy the constraints of equation [3] is ~ 23 
M Q . Assuming the inclination of the system is close to 90 deg, the mass of the black hole 
is most plausibly in the range 23 — 34 M Q . However, if the inclination is much lower than 
expected (close to 45 deg) the black hole mass can be as high as ~ 57 M . 

If the mass of the the compact object in IC 10 X-l is indeed >23 M , it is the most 
massive known stellar black hole with a dynamically a determined mass (the record is 
curently held by GRS 1915+105 with a mass 10-18 M fjllemillard fc McClintockl bood ).) 
Mos t models predict that the remnants of a supernova explosion have masses below 20 
M ( lFryer fc Kalogerall200ll ). It is possible that the bl ack hole in IC 10 X-l has increased 
it's mass by accreting material from [MAC92] 17-A (iPodsiadlowski et al.l 120031 ) The low 
metallicity of IC 10 ma y also be conducive to the formation of higher mass black holes 
( IFryer fc Kalogerall200ll ). Detailed theoretical models need to be developed to understand 
the evolutionary history of IC 10 X-l. 



3.3. Is the Model of IC 10 X-l Self Consistent? 

The results from this and previous studies of the IC 10 X-l system suggest that it is 
a detached binary composed of a massive (~ 35 M ) Wolf-Rayet star with a 23-34 M 
black hole companion. X-rays are produced as material from the wind of [MAC92] 17-A is 
accreted via an accretion disk onto the black hole. It is important to test the self-consistency 
of this model and show that such a system is indeed capable of forming an accretion disk 
and producing X-rays. 

A black hole in a detached binary system can form an accretion disk (and hen ce be a 



strong X-ray source) if the following condition is met e.g. tjErgma fc Yungelsonl 119981 ): 



P < A.8(M BH /M Q )v^ 00 5 2 (4) 

where Mbh is the mass of the black hole, v the velocity of the wind impacting the compact 
object in units of 1000 km s _1 , P is the period in hours, and 5 a dimensionless parameter of 
order unity. If the mass of the black hole is 24 — 34 M we find that this condition is satisfied 
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if v < 1400 — 1500 km s 1 . The terminal wind velocity of [MAC92] 17- A was estimated by 
Clark fc Crowtherl (120041 ) to be = 1750 km s" 1 . The velocity of the wind impacting the 



black hole will be less than the terminal velocity; the up per limit of 1400 — 150 km s 1 is 



consistent with the wind structure of a Wolf-Rayet star (ICarpano et al. Il2007bl ). 



The energy released by material falling onto a compact object via an accretion disk 
is given by L = rjMc 2 , (IShakura fc Syunyaevl Il973l ). where M is the accretion rate and rj 
the efficiency If we assume r] = 0.1, the mass accretion rate required to sustain an X-ray 
luminosity of L x ~ 2 x 10 38 ergs" 1 i s ~ 3. 5 x 10" 8 M yr _1 . The mass loss rate of [MAC92] 
17- A derived by Oark fc Crowtherl (120041 ) is ~ 10~ 5 M Q yr _1 . We therefore conclude that 



direct wind fed accretion is consistent with the observed X-ray luminosity. 



3.4. IC 10 X-l and the Star Formation History of IC 10 



Both the mass of the donor and the mass of the compact object in IC 10 X-l indicate 
that this system is less than 10 Myr old. This is consistent with the suggestion th at there 
have been pockets of star formation ~ 3 — 4 Myr ago which may still be on-going (IHunter 
200 11). There are several dozen WR stars in IC 10, also indicating recent star formation 
( IMassey fc Holmesll2002l ; ICrowther et al.ll2003l ). The fact that the donor is a He- core WR star 
i ndica tes that its initial mass was much higher (probably in excess of 60M g lMeynet fc Maeder 
( 120051 )) in which case the system could have formed in one of the recent star formation events. 
The presence of a shell of ionized gas around the system, probably related to supernova 
activity, indicates that there is still on-going star-formation in that region. The f ormation of 
the very massive progenitors of this binary requires a very flat IMF. iHunterl (120011 ). finds that 
such a flat IMF is consistent with a star-formation scenario of widespread activity ~ 40 Myr 
ago and more recent (~ 3 — 4 Myr) localized events. 



4. Summary and Conclusions 



We have used SWIFT and Chandra to determine the orbital period of the bright Wolf- 
Rayet X-ray binary IC 10 X-l to be 34.4 hours. Us i ng th is value of the period, and the 
mass of the donor star derived by IClark fc Crowtherl (120041 ) from optical spectroscopy, we 
use Kepler's laws to show that the donor star almost certainly does not fill its Roche lobe. 
We cannot derive the mass function of the binary system or mass of the compact companion 
from the period alone. However, re-examination of published optical data reveals a shift in 
the centroid of the Hell 4686A line most likely due to orbital motion of the donor star. If 
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this is the case, the mass function f(M) ~ 7.8 M and the compact object is a black hole. 
Assuming that the mass of the donor star is 17 — 35M M , we calculate that the black hole 
mass is ~ 24 — 34 M Q . If confirmed, this makes IC 10 X-l the most massive black hole binary 
known. We show that it is possible for an accretion disk to form and that the observed mass 
loss rate can power the X-ray luminosity. 
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Table I. Black Hole Mass (in M ) as a Function of Inclination and Donor Mass 





Donor Mass (M ) 




17 25 35 


Inclination (deg) 




90 


23 28 34 


60 


29 35 41 


45 


43 49 57 
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Fig. 1. — SWIFT lightcurve of IC 10 X-l in 100 second time bins. Overplotted in red is 
the Chandra lightcurve in 1000 second bins. The Chandra count rate is reduced by a factor 
of 5 to approximate the lower effective area of SWIFT. The Chandra data is also shifted 
forward in time by 13 x p where p is the Lomb-Scargle period. Times are relative to the 
start of S WIFT observations, 2006 21 November at 05:11:05 UTC. The blue curve is the 
period derived by the Lomb-Scargle periodogram. 




O 50 100 150 200 250 
Period [hours) 



Fig. 2. — The power spectrum of IC 10 X-l. The horizontal dotted lines correspond to the 
false- alarm probability corresponding to 1, 2, and 3-sigma significance. 




Fig. 3.— The Hell 4686 Aline in [MAC92] 17- A. Four individual spect ra are shown. Zero 
velocity is shown as the dashed line. There is clear evidence that the line centroid has shifted. 



